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ABSTRACT Positional cloning in maize (Zea mays) requires development of markers in the region of in- 
terest. We found that primers designed to amplify annotated insertion-deletion polymorphisms of seven 
base pairs or greater between B73 and Mo1 7 produce polymorphic markers at a 97% frequency with 49% of 
the products showing co-dominant fragment length polymorphisms. When the same polymorphisms are 
used to develop markers for B73 and W22 or Mo17 and W22 mapping populations, 22% and 31% of 
markers are co-dominant, respectively. There are 38,223 Indel polymorphisms that can be converted 
to markers providing high-density coverage throughout the maize genome. This strategy significantly 
increases the efficiency of marker development for fine-mapping in maize. 
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Maize whole genome sequencing has greatly simplified positional 
cloning of mutant loci and quantitative traits. Maize has very high 
levels of nucleotide sequence diversity (Chia et al. 2012; Jiao et al. 

2012) , which has been leveraged to develop multiple molecular 
marker platforms (Sharopova et al. 2002; Fu et al. 2006; Gore et al. 
2009; Liu et al. 2010; Frascaroli et al. 2013; Qu and Liu 2013; Xu et al. 

2013) . These marker sets enable mapping of phenotypic loci, but the 
primary challenges to positional cloning are efficient identification of 
recombinant chromosomes and rapid development of molecular 
markers in the region of interest. A typical positional cloning project 
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begins by using one of the developed and low-cost mapping platforms 
that rely on single nucleotide polymorphisms (SNPs) to map the locus 
(Jander et al. 2002; Liu et al. 2010). A large mapping population is 
then screened to identify recombinant individuals within the region. 
Finally, high-density markers are needed to narrow the area to a few 
candidate genes. 

Flexible and scalable marker systems utilizing simple sequence 
repeat (SSR) or insertion-deletion (Indel) polymorphisms are a com- 
mon choice for fine-mapping a single locus (Phillips et al. 2011; 
Slewinski et al. 2012). SSR tracks can be identified in the B73 genome 
sequence, but these require experimental testing to determine if in- 
dividual markers will be useful for different mapping applications 
(Martin et al. 2010). Indel markers can be identified at low frequency 
by amplifying PCR products spanning introns or 3' UTR of annotated 
genes (Fu et al. 2006). 

The Mo 17 inbred was sequenced by the Department of Energy 
Joint Genome Institute (JGI) with SNP and Indel polymorphisms 
annotated on the B73 genome (Sen et al. 2009). Markers designed 
from these annotations should be polymorphic for B73/Mol7 map- 
ping populations and are expected to have a high frequency of poly- 
morphism when either B73 or Mol7 is used as a mapping parent. 
Here, we show that PCR primers designed to amplify annotated B73/ 
Mol7 Indel polymorphisms produce a high frequency of co-dominant 
and dominant molecular markers that can be scored using agarose gel 
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Figure 1 Examples of Indel markers designed from the Mo17 poly- 
morphism track of MaizeGDB. (A) Marker uflDPI -256.9 shows the most 
common co-dominant marker class, with B73 and W22 sharing the same 
size allele and each having a co-dominant polymorphism with Mo17. The 
B73 product was designed to amplify a 1 54-bp product, with Mo1 7 having 
a 17-bp deletion. (B) The uflD6-1 33.83 marker shows an Indel that 
amplified three different co-dominant alleles with a 255-bp B73 product 
and a 19-bp Mo17 deletion. (C) The uflDP7-16.1 marker shows a co- 
dominant polymorphism between B73 and Mo17, whereas W22 shows 
PAV polymorphisms with both B73 and Mo17. The B73 product was 
designed to amplify a 1 43-bp product with a 1 3-bp Mo1 7 deletion. (D) The 
uflDPI -257.3 marker amplifies a 168-bp product from B73 and W22 with 
a 14-bp insertion in Mo17. The B73 and W22 products are dominant over 
the Mo17 product. (E) The uflDP4-20.905 marker amplifies a 201-bp 
product from B73 with an 11 -bp deletion in Mo17 and W22. The B73 allele 
is dominant. (F) The uflDP4-25.42 primers amplify a 192-bp product from 
all three inbred lines despite a 24-bp deletion predicted for Mo17 based 
on the polymorphism track. 

electrophoresis. This marker development strategy enables rapid 
fine-mapping of mutants in any laboratory equipped for basic molec- 
ular biology techniques. 



MATERIALS AND METHODS 

MaizeGDB (www.maizegdb.org) provided the complete Mo 17 poly- 
morphism track based on the JGI sequence, which contains more than 
4.5 million SNP and Indel polymorphisms. We selected the subset of 
38,223 Indel polymorphisms with at least a 7-bp difference between 
the B73 and Mol7 alleles to design PCR primers for fine-mapping 
(Supporting Information, Table SI). Primer pairs were designed to 
amplify PCR products spanning 330 annotated Indels with expected 
B73 PCR products ranging from 80 to 300 bp in length. All primers 
selected were specific to the target B73 locus based on BLASTn 
searches of the maize genome assembly release 5b.60 and were tested 
using inbred DNA from B73 and Mo 17. Subsets of primer pairs 
were tested for amplification of W22 inbred as well as F x or mixed 
DNA from B73/Mol7, B73/W22, or Mol7/W22 to determine if 
the markers were co-dominant (Figure 1). PCR and 4% agarose gel 
electrophoresis were completed as described by Martin et al. 
(2010). 

RESULTS AND DISCUSSION 

We tested 330 markers for mapping 15 genetic loci to sample nearly 
334 Mbp of the genome (Table S2). A B73 product amplified from all 
primer pairs tested, and a Mol7 allele amplified for 240 of the 330 
markers. The PCR primer design strategy did not account for nearby 
Mol7 SNPs and short Indels. It is possible that Mol7 alleles failed to 
amplify in 90 loci due to nucleotide divergence between the B73 
primer sequences and the Mol7 loci targeted. For the 240 markers 
where both B73 and Mol7 alleles amplified, we evaluated whether the 
expected size differences were observed based on the annotated in- 
sertion or deletion (Figure 1). Thirty markers amplified unexpected 
products from Mol7. In 12 cases, deletion polymorphisms amplified 
a larger Mo 17 product or insertion polymorphisms amplified a smaller 
Mol7 product. For the remaining 18 markers, the Mol7 and B73 PCR 
products were the same size (Figure IF). These data indicate a low 
level of disagreement (10% of Indel annotations) between the genome 
annotation and experimental evidence. 

We tested 259 markers with B73/Mol7 DNA (Table 1). Of these, 
49% amplified both B73 and Mo 17 alleles as co-dominant fragment 
length polymorphisms (Figure 1, A-C). The small size differences 
between the B73 and Mol7 alleles resulted in frequent formation of 
slow migrating, heteroduplex bands in the B73/Mol7 PCR products. 
An additional 48% of the markers amplified only one allele in the B73/ 
Mol7 DNA and were scored as dominant or presence-absence var- 
iants (PAV) (Figure 1, D and E). In all but one case, the B73 allele 
preferentially amplified over the Mo 17 allele, which is consistent with 
the primers being designed from the B73 reference genome. Only 3% 
of markers were not polymorphic (Figure IF). 

Fu et al. (2006) reported a strategy to generate Indel markers 
between B73 and Mo 17 by amplifying intron or 3' UTR sequences. 
Slightly more than 7% of nearly 12,000 PCR products showed PAV 
polymorphism, and an additional 3.5% of the products showed dis- 
tinct fragment length polymorphisms. This latter set of length poly- 
morphisms was not tested for co-dominance in Fj or mixed DNA. 
A subset of these markers would be expected to show dominant am- 
plification similar to that shown in Figure 1, D and E. Thus, designing 



Table 1 Summary of Indel markers tested using B73/Mo17, B73/W22, and Mo17/W22 DNA 



Inbred Pair 


No. Markers Tested 


Co-Dominant 


PAV 


Dominant 


Not Polymorphic 


% Co-Dominant 


B73/Mo17 


259 


127 


80 


44 


8 


49 


B73/W22 


201 


44 


36 


16 


105 


22 


Mo17/W22 


201 


63 


48 


20 


70 


31 
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markers from the Mol7 polymorphism track yields at least 14-fold 
higher frequency of co-dominant markers than the Fu et al. (2006) 
strategy. 

Many public mutagenesis resources utilized the W22 inbred 
background for mutagenesis (Cowperthwaite et al. 2002; Till et al. 
2004; Kolkman et al. 2005; McCarty et al. 2005; Ahern et al. 2009). 
To determine how readily B73/Mol7 Indel polymorphisms can be 
applied to fine-mapping with a W22 genetic background, we screened 
201 markers for co-dominant polymorphisms in both the B73/W22 
and Mol7/W22 inbred pairs (Table 1). The W22 allele was frequently 
the same length as the B73 allele with 52% of markers amplifying the 
same size product (Figure 1A), whereas 35% of the markers amplified 
the same size allele from both W22 and Mo 17 DNA. Polymorphic 
markers for these two sets of inbred parents were divided equally 
between co-dominant and dominant/PAV polymorphisms (Table 
1). Overall, 86% of the 201 Indel markers tested for B73AV22 and 
Mol7/W22 were polymorphic with at least one pair. Compared to 
SSR markers, only 43% of 4083 randomly selected SSR tracks were 
found to be polymorphic among 11 diverse inbred lines (Sharopova 
et al. 2002). More recentiy, next-generation sequencing of hundreds of 
inbred lines has identified thousands of SSRs with length polymor- 
phisms predicted at a similar frequency to that observed by Sharopova 
et al. (2002) (Qu and Liu 2013; Xu et al. 2013). To identify SSR loci 
with as high a frequency of polymorphic markers as found in the B73/ 
Mo 17 Indel annotation, specific comparisons between mapping pop- 
ulation parents would be needed. 

Based on our review of current public resources, we believe Table 
SI provides polymorphisms with the highest likelihood for successful 
development of new co-dominant amplified fragment length markers 
as long as either the B73 or the Mo 17 inbred is a parent for the 
mapping population of interest. The average distance between potential 
markers in Table SI is 53.7 kb, with a median distance of 10 kb 
between polymorphisms. Larger distances between polymorphisms co- 
incide with centromere and heterochromatic regions of the chromo- 
somes. Thus, polymorphism density is highest in gene-rich regions, 
and it is expected that the polymorphisms in Table SI should provide 
sufficient coverage to fine-map any trait to a few candidate genes. 

No trends were found for B73/Mol7 Indel marker co-dominant 
polymorphism rates based on the size of the Indel or the size of the 
expected B73 product. Heteroduplex products allow heterozygous 
recombinant individuals to be readily scored on agarose gels even if 
the Indel has a small size difference. However, a smaller (<150 bp) 
PCR product will make smaller size differences (<10 bp Indel) be- 
tween contrasting homozygous alleles easier to resolve in agarose gels. 
Key design concerns for successful marker development are ensuring 
that the PCR primer sequences are specific for the target Indel locus 
and that the expected PCR product is within the range of resolution 
for a 4% agarose gel. To further improve the successful development 
of co-dominant markers, primer design could also take into account 
linked Mo 17 SNP and short Indel polymorphisms. Our current fine- 
mapping strategy is to select five evenly spaced polymorphisms in the 
mapping interval to identify one to three co-dominant markers that 
refine the interval. We then repeat the process with evenly spaced 
polymorphisms within refined intervals until a suitably small region 
is identified for sequencing candidate genes. 
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